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Abstract 

We present some formulae for the moments of inertia of ho- 
mogeneous solids of revolution in terms of the functions that 
generate the solids. The development of these expressions ex- 
ploits the cylindrical symmetry of these objects, and avoids 
the expHcit use of multiple integration, providing an easy and 
pedagogical approach. The expUcit use of the functions that 
generate the solid gives the possibility of writing the moment 
of inertia as a functional, which in turn aUows us to utiHze 
the calculus of variations to obtain a new insight into some 
properties of this fundamental quantity. In particular, min- 
imization of moments of inertia under certain restrictions is 
possible by using variational methods. 

Keywords: Moment of inertia, variational methods, 
solids of revolution. 

PACS: 45.40.-F, 46.05.th, 02.30.Wd 

The moment of inertia (MI) is a very important concept in 
Physics and Engineering 0. In this paper, we present simple 
formulae to obtain the Mi's of homogeneous solids of revolution. 
The expressions presented here are written in terms of the func- 
tions that generate the solid and only require simple integration. 
Finally, we show that minimization of the moment of inertia un- 
der certain restrictions is possible by employing the calculus of 
variations. 

1 MI of solids of revolution gener- 
ated around the X— axis 

Figure ^ shows a function that generates a solid of revo- 
lution around the X— axis (we shall call it the "generating 
function" henceforth). The narrow rectangle of height f (x) 
and width dx generates a thin disk of height dx and radius 
f (x). We shall calculate the moment of inertia of the solid 
of revolution generated by / (x) with respect to the axis of 
symmetry (X— axis). 

We know from the literature [D , that the Ml of a thin disk 
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Figure 1: Solid of revolution generated from the X~axis. (a) 
With one generating function f (x). (b) With two generating 
functions fi (x) , /2 (x). 

with respect to the X— axis is given by (1/2) MR^ where M 
is the mass of the disk and R gives its radius. Thus, the 
differential Ml for our thin disk reads 

I 



dix = - [dM] f {xY . 

The differential of mass is given by 

dM = p dV = pnf (a;)^ dx, 



(1) 



(2) 



where p denotes the density of the solid, and will be assumed 
constant throughout the document. Substituting Eq. Q, 
into Eq. 0J and integrating we get 

=/ 
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f{x) dx. 



(3) 



Eq. (PI gives the MI of any solid of revolution with respect 
to the axis of symmetry. Now let us calculate the MI with 
respect to the y— axis. To do it, we first estimate the MI of 
the thin shaded disk around the axis Y' shown in Fig. QJi. It 
is well known that for an axis passing through the diameter 
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of the disk, the MI reads ly = (1/4) Mi?^ [T], in our case 
liave 1 „ , .9 dix 



dly 



{dM) f {xY 



(4) 



4 ' '■' ' ' 2 ' 
to calculate the differential MI with respect to the F— axis 
we use the parallel axis theorem obtaining 



dly = dly + x^dM 



dl 



X 



x^dM. 



Replacing Eq. (0) into Eq. jSj) and integrating we find 



lY = -^+7rp[ x^f{xfdx 



(5) 



(6) 



This expression provides the perpendicular MI of the solid * . 
The perpendicular Mi's are not usually calculated in common 
textbooks. Notwithstanding, they are important in many 
physical problems. For instance, some solids of revolution 
acting as a physical pendulum require this perpendicular MI. 
Eqs. (0 EJ show that Ix and /y can be calculated with 
a simple integral based on the generating function of the 
solid. This simplification comes from the cylindrical or axial 
symmetry that solids of revolution exhibit. Finally, axial 
symmetry also tells us that ly = Iz- 

By examining Eq. ISJ we find that dly' = dlx/"^, it comes 
from the combination of the perpendicular axis theorem and 
the axial symmetry applied to the thin disks. Further, the 
term x'^dM in Eq. lO comes from the parallel axis theorem. 
These facts give a geometrical interpretation of Eq. © , the 
first term on the right-hand side of such an equation comes 
from the combination of the perpendicular axis theorem with 
the axial symmetry applied to the element of volume, while 
the second term comes from the parallel axis theorem. 

On the other hand, if we are interested in solids of revo- 
lution generated by two functions /i (x) and /2 {x) as Fig. 
^ displays, we only have to substract the contribution of a 
solid generated by /i (x) from the figure generated by /2 (x) 
getting 
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dx. 



(7) 



h{xY-h{xy dx. (8) 



We have assumed that /2 (x) > /i (x) for x £ [xq , x/] . Once 
again, Eqs. Q |SJ) only require simple integration and the 
knowledge of the generating functions. It is worth pointing 
out that all the development given in this section, still holds 
if the density depends on the x variable only i.e. on the 
"height" of the solid, except that in the expressions (0 El 13 
inj, the density should be inside the integrals. It comes from 
the fact that when p = p{x), each one of the thin disks used 
in the demonstration is still homogeneous. 

*For "the perpendicular MI" wc mean the MI with respect to an axis 
perpendicular to the axis of symmetry. There are an infinite number of 
such axes but all these Mi's are related by the parallel axis theorem. 




X 



Figure 2: Frustum of a right circular cone with a conical well. 
The shadowed surface is the one that generates the solid, and 
fi (x), /2 (x) are the generating functions. 

Example 1. Mi's for a truncated cone with a conical well 
(see Fig. The generating functions read 



/i(^-) = 



Ril 



if X e [0, h] 

if X £ {h, H] 



X + a2 



(9) 



where all the dimensions involved are displayed in Fig. |5| Substi- 
tuting Eqs. into Eqs. (QISJ we get 
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Ix = ^[H [ai + a2 + aia2 + aia2 + aia2) — R h] (10) 



+ 



R^ 
6 



(11) 



It is more usual to give the radius of gyration (RG) instead 
of the MI. For this we calculate the mass of the solid, whose 
expression in terms of the generating functions is well known from 
the literature 



M 



np 



[/2(a 



Y-fi{xf] dx 



(12) 



from Eq. iHH . we get M and the RG's become 

3 {H (af + ai+ aiai -I- a?a2 + a?ai) 



K'x 



R'h} 



10 [H (aiaa + + a%) - K^h 
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[H (aiaa + a'f + a^) - R-^h] 



(13) 



By setting R — (and/or h = 0) we find the RG's for the 
truncated cone. With R = and ai = 0, we get the RG's of a 
cone for which the axes Y and Z pass through its base. Setting 
R = and 02 = 0, we find the RG's of a cone but with the axes 
Y and Z passing through its vertex. Finally, by setting R = 0, 
and ffli =02; we obtain the RG's for a cylinder. In many cases 
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of interest, we need to calculate the Mi's for axes Xc, Yc and 
Zc passing through the center of mass (CM), these Mi's can be 
calculated by finding the position of the CM with respect to the 
original coordinate axes, and using the parallel axis theorem. The 
position of the CM can be easily found to be {xcm, 0, 0) with 



XCM = 



/;/ X [h {xf - h i^f] dx 

[h i^r - h (^)l dx ■ 



(14) 



Applying Eq. (I14II the position of the CM for the truncated cone 
with a conical well reads 



XCM ~ 



[(2aia2 + 3a? + al) - R^h" 



4 [H (aia2 + a? + a^) 



Gathering Eqs. ^ we find 
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Kt 
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(15) 



(16) 



2 MI of solids of revolution gener- 
ated around the F— axis 
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Figure 3: Solid of revolution generated from the Y~axis. 

A generating function can be used to form a solid from the 
y— axis as Fig. |3] indicates. In order to calculate the MI of this 
figure with respect to the F— axis, we should calculate the differ- 
ential of mass, dM — pdV, corresponding to the cylindrical shell 
shown in Fig. |3] 



dM = pTTf{x)[{x + dx) — X ] = 2npxf{x) dx, 



(17) 



where differentials of second order are neglected. It is clear that 
dlv = x^dM, and by integrating we find 



Iy = 2tvp 



X f {x) dx. 



(18) 



It gives the MI of the solid with respect to the axis of symmetry 
(y— axis). In order to calculate the MI around the X-axis, we first 
estimate the perpendicular MI (with respect to an axis passing 
through the CM) of a homogeneous cylindrical shell with inner 
radius ai, outer radius a2, and height h. It can be calculated by 
replacing xq = —h/2, Xf = h/2, with f2{x) — a2 and fi{x) = ai 
in Eq. JSJ obtaining ^ 

tit is important to take into account that in Sec. Q the axis of 
symmetry is the X-axis, while in section ||2J the axis of symmetry 
is the y— axis. 



iXcM = — l«2 ~ oil) " + 



TTp I aa 



a?) 



For our particular cylindrical shell we have ai = x, a2 = x + dx, 
h = f (x); from which Ixcm becomes differential. Neglecting 
differentials of second order we find 

\3 

dx. (20) 



3 TV pX f (x) 

dIx,CM = T^px f (x) dx H • 



6 

Since in general each infinitesimal cylindrical shell has a different 
center of mass, we cannot integrate this result directly to obtain 
IxcM- Instead, we shall use the parallel axis theorem to find dix, 
i.e. the MI of the cylindrical shell with respect to the X— axis. 
From Eqs. 1171 1201 and using the parallel axis theorem, we get 



dIx =dIx,CM + dM 



2 

= npx^f (x) dx + -Tvpx f (x)^ dx. 



Integrating in x and taking into account Eq. H18|l gives 

X f (x) dx. 



(21) 



(22) 



If the solid is generated by two functions f2{x) and /i(a::) we 
can make a substraction like in the previous section, and Eqs. 1181 
I22I I become 



Iy = 2'Kp 



J xn 



[/a (x) - fi {^)] dx. 



Y + ^l ^ x[f2{xf-f^ixf] dx. 



(23) 

(24) 



When the figure is generated around the y— axis, we should as- 
sume that xo > 0; such that all points in the generating sur- 
face always have non-negative x coordinates. Instead, we might 
allow fi (x), f2 (x) to be negative though still demanding that 
/i {^) 5- /a (x) in the whole interval of x. Once again, the cylin- 
drical symmetry indicates that Ix = Iz. 

As in the previous section, the expressions for the mass and the 
center of mass for figures generated around the y— axis can also 
be derived easily 



M =27Vp / X [/2 (x) - /i {x)] dx 



ycM 



J X [/2 {xf - fi {xf] dx 
' 2jx[f2{x)-fi{x)] dx • 



(25) 



these expressions are important to calculate RG's and Mi's around 
axes passing through the CM of the figure. 

Equations 123II24II are especially useful in the case in which the 
generating functions /i (x), /a (a::) do not admit inverses, because 
in such a case we cannot find the corresponding inverse functions 
gi (x), g2 {x) to generate the same figure by rotating around the 
X— axis. This is the case in the following example 

Example 2. Calculate the Mi's of a solid formed by rotating 
the function f {x) — h + A sin {nnx / R) , around the Y — axis (see 
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Figure 4: Solid of revolution created by rotating the generat- 
ing function f (x) = /i+Asin(^2^) around the Y— axis. From 
the picture it is clear that f{x) does not admit an inverse. 

Fig- IDi where the function is defined in the interval x G [0,ii], 
and n is a positive integer. We demand /i > if n > 1; besides, 
if n = 1 and \A\ > h we demand ^ > 0. These requirements 
assure that f (x) > for all x € [0,i?]. Replacing f{x) into Eqs. 



1181 \'2'2l and calculating the mass with Eq. 
RG's 

i?2 ^n7r/i + 4A(6n-^7r-2 - !)(-!)' 



I25II we obtain the 



+ 



n7r/i + 2A(-l)"+i 
3n-Kh[2h'^ + SA^] + 4A{- 



(26) 



2 ' 18[n7r/i + 2A(-l)"+i] 

Observe that / (x) does not have an inverse. Hence, we cannot 
generate the same object by constructing an equivalent function 
to be rotated around the X— axis. This figure could for instance, 
simulate a solid of revolution with a rugged surface on the top, or 
the furrows formed in a piece of material that has been machined 
in a lathe. 




Figure 5: A bell formed by two Gaussian distributions rotat- 
ing around the Y—axis, the bell tolls around an axis perpen- 
dicular to the axis of symmetry, that passes through the point 
B 

Example 3. MPs for a Gaussian Bell. Let us consider a hollow 
bell, which can be reasonably described by a couple of Gaussian 
distributions (see Fig. |^. 



f^{x) = Ae- 



f2{x) = Be 



(27) 



where a, (3, A, B are positive parameters (0 < ^4 < B, and a > /3 
) that allow us to model the Bell, B is the height of the Bell, 
B — A is its thickness on the top, 1/q and 1//3 are the decays 
that simulate the profile. For the sake of simplicity, we integrate 
fl (x) , /2 (2;) in the interval (0,oo), however integration in a 
finite interval does not change the results significantly if the bell 
is wide enough. The MPs are obtained from 1231 1241 



Iy ~ vrp 



B_ 



= -f 



7rp 


■ B 
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B' 


A'- 
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The mass and the center of mass position read 



M : 



np 



ycM 



aB^ - PA^ 



aB - f3A 



xcM = . (29) 



When the bell tolls, it rotates around an axis perpendicular 
to the axis of symmetry that passes the top of the bell. Thus, 
this is a real situation in which the perpendicular MI is required. 
On the other hand, owing to the cylindrical symmetry, we can 
calculate this MI by taking an axis parallel to the X— axis that 
passes through the top of the bell, which corresponds to y = B. 
By using the parallel axis theorem it can be shown that 

Ix,B =Ix+ MB{B - 2Ycm) 



18a2/32 



[q^B(9 + llB^/3) 



+ l3'^A{'SABa - 2A^a - 18B^a - 9)]. 



(30) 



This example shows another advantage of working with a solid 
generated around the F— axis. If you try to do the same problem 
by employing Eqs. Q EJ you need the inverse of the gaussian 
distribution and the integration is more complex. Indeed, in many 
cases integration involving the inverse of certain function could be 
much harder than integrations involving the function itself, even if 
the inverse exists. This is another justification to develop formulae 
for solids generated around the axis. 

3 Applications utilizing the calculus 
of variations 

In all the equations shown in this paper, the Mi's can be seen as 
functionals of some generating functions [5]. For simplicity, we 
take a solid of revolution generated around the X— axis with only 
one generating function. From Eqs. (0 El we see that the Mi's 
are functionals of / (x), so that 



Ixlf] 
Ivif] 



np 
~2 

Ix[f] 



f{x'Ydx', (31) 
- -Kp j x'^ f (s') dx' . (32) 



Then, we can use the methods of the calculus of variations (CV) 
in order to optimize the MI. To figure out possible applications, 
imagine that we should build up a figure such that under certain 
restrictions (that depend on the details of the design) we require 
a minimum of energy to set the solid at certain angular velocity 
starting from rest. Thus, the optimal design requires the moment 
of inertia around the axis of rotation to be a minimum. 

As a specific example, suppose that we have a certain amount of 
material and we wish to make up a solid of revolution of a fixed 



(28) 



tThe reader not familiarized with the methods of the CV, could 
skip this section without sacrificing any understanding of the rest 
of the content. Interested readers can look up in the extensive 
bibliography concerning this topic, e.g. Ref. [3- 



Moments of inertia for solids of revolution and variational methods 



5 



length with it, such that its MI around a certain axis becomes 
a minimum. To do it, let us consider a fixed interval [a;o,a;/] of 
length L, to generate a solid of revolution of mass M and constant 
density p (see Fig. ^). Let us find the function / (x), such that 
Ix or ly becomes a minimum. Since the mass is kept constant, 
we use it as the fundamental constraint 

M = TTp f (a-') dx' = constant. (33) 

JXQ 

In order to minimize Ix we should minimize the functional 



) 










m 


L/2\ ' 



-R/fi R/fi \ y 

Figure 6: (a) Elliptical function that generates the solids of 
revolution that minimize ly ■ The shaded region is the one 
that generates the solid, (b) Truncated spheroid obtained 
when the shaded region is revolved. 



[/] = r 9 if,^') dx' = Ix [f]-Xnp r f {x'f dx' (34) 

J xn ^ xn 



where A is the Lagrange multiplicator associated with the con- 
straint 13311 . In order to minimize Gx [f], we should use the 
Euler-Lagrange equation [3| 



SGx[f] _ dg{f,x) d dg{f,x) 



Sfix) Of dxdidf/dx) 







obtaining 



SGxlf] 
Sf{x) 



27vpf {xf ^ 27v\pf (x) = 



whose non-trivial solution is given by 

f(x) = ^/\ = R 



(35) 

(36) 
(37) 



Analysing the second variational derivative we realize that this 
condition corresponds to a minimum. Hence, Ix becomes mini- 
mum under the above assumptions for a cylinder of radius \/A, 
such radius can be obtained from the condition 13311 . yielding 
= M/ TvpL and Ix becomes 



- iii!. 

J- X .cylinder "7 7" 

2 npL 



(38) 



Now, we look for a function that minimizes the MI of the solid of 
revolution around an axis perpendicular to the axis of symmetry. 
From Eqs. (ElOSJ, 

we see that the functional to minimize is 

Gr[.f] = ^^ + .P 



x'^ f {x'^'^ dx —X-Kp I f (^x')^ dx' , 

(39) 

setting the variation of Gy [/] with respect to / (a;) we get 



/ {xf =2{\-x^)=R^- 2x^ , (40) 

where we have written 2A = R^. By taking Xq = — iv/2, Xf = 
L/2, the function obtained is an ellipse centered at the origin with 
semimajor axis R along the axis, semiminor axis R/^/2 along 
the X— axis, and with eccentricity e = 1/ ^/2. When it is revolved 
we get an ellipsoid of revolution (spheroid); such spheroid is the 
solid of revolution that minimizes the MI with respect to an axis 
perpendicular to the axis of revolution. From condition 1331 we 
find 



M 

TTpL 
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(41) 



In the most general case, the spheroid generated in this way 
is truncated, as shown in Fig. |H] and the condition R > L/\/~2 
should be fulfilled for / (x) to be real. The spheroid is complete 
when R = and the mass obtained in this case is the mini- 

mum one for the spheroid to fill up the interval [—L/2, -L/2], this 
minimum mass is given by 



from (El, m, (EIJ and ^ we find 

M^inL^ [bp + - 1] 



M 



'-Y, spheroid 



60 



P ■: 



(42) 



(43) 



Assuming that the densities and masses of the spheroid and the 
cylinder coincide, we estimate the quotients 

lY,cylinder _ (^A* + — l) 



lY,Bpheroid 5/i {p + 1) 



^ X ,cylinder 
I X ,spheroid 



= 1 



1 -2 



< 1, 



< 1 



(44) 



Eqs. I^J show that ly^sph < lY,cyi while Ix,cyi < Ix,sph- In 
both cases if M >> Mmin the Mi's of the spheroid and the cylin- 
der coincide, it is because the truncated spheroid approaches the 
form of a cylinder when the amount of mass to be distributed in 
the interval of length L is increased. 

On the other hand, in many applications what really matters 
are the Mi's around axes passing through the CM. In the case 
of homogeneous solids of revolution the axis that generates the 
solid passes through the CM, but this is not necessarily the case 
for an axis perpendicular to the former. If we are interested in 
minimizing ly^,, i.e. the MI with respect to an axis parallel to 
Y and passing through the CM, we should write the expression 
for Jy^ by using the parallel axis theorem and by combining Eqs. 



Ix[f] 



np 



TTp 



Qfix'Ydx' 



x'^ f (a;') ^ dx' 



x f{x')^ dx' 



(45) 



thus, the functional to be minimized is 
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Gyc [f] = Iyc [/] - A^p r f {xf dx' , (46) 

J Xg 

after some algebra, we arrive to the following minimizing function 
f{xf^R^-2{x-xcMf , (47) 

where we have written 2 A = R^. It corresponds to a spheroid 
(truncated, in general) centered at the point (a;cA/,0,0) as ex- 
pected, showing the consistency of the method. 

Finally, it is worth remarking that the techniques of the CV 
shown here can be extrapolated to more complex situations, as 
long as we are able to see the Mi's as functionals of certain gen- 
erating functions. The situations shown here are simple for ped- 
agogical reasons, but they open a window for other applications 
with other constraints^, for which the minimization cannot be 
done by intuition. For instance, if our constraint consists of keep- 
ing the surface constant, the solutions are not easy to guess, and 
should be developed by variational methods. 



Engineering Mechanics Statics, 7th Ed. (Prentice-Hall Inc., 
New York, 1995). 

[2] For additional examples concerning the formulae 
shown here, the interested reader is referred to the 
website: xxx.lanl.gov /list /physics / 0404 
Preprint number: physics/0404005J 

[3] G. Arfken. Mathematical methods for physicists. Second Ed. 
(Academic Press, 1970) Chap. 17. 



4 Conclusions 

We found some formulae that provides a direct and simple way 
to calculate moments of inertia for homogeneous solids of revolu- 
tion. The approach is easy and pedagogical to be used in basic 
courses of Physics and Engineering, because multiple integration 
is avoided in the demonstration by taking advantage of the cylin- 
drical symmetry that solids of revolution possess. It is worth 
emphasizing that perpendicular moments of inertia are not usu- 
ally reported in the literature although they are also important in 
many physical applications; however, Eqs. llSl l241 show that they 
can also be evaluated by simple integration. In addition, Eqs. ((3 
IH1ES1I2U show that we do not have to worry about the partitions 
usually required when multiple integration is done, we only have 
to know the generating functions of the solid. 

On the other hand, it deserves to point out that the formulae 
shown here can be seen as functionals of the generating functions 
of the solids. This fact permits the use of the calculus of varia- 
tions to explore the properties of moments of inertia. In particu- 
lar, minimization of moments of inertia under certain restrictions 
is possible by using variational methods. It could be useful in 
applied Physics and Engineering. 

The authors acknowledge to Dr. Hector Miinera for revising 
the manuscript. 
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